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Abstract 


Approximate formulas are derived for the intensity of diatomic molecular elec- 
tronic transitions. The formulas correspond to various models of the spectrum, 
ranging from a monochromatic approximation to others in which more details of 
the vibrational and rotational structure are taken into account. Some of these 
methods are well known, while others appear to be new. Choice of the appropriate 
method in a given circumstance depends upon the spectral shape, knowledge of 
the molecular parameters, and the desired accuracy of the calculation. Each model 
is applied to the CN violet system to illustrate its accuracy. 
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Approximate Formulas for the Intensity of Electronic 
Transitions in Diatomic Molecules 


I. Introduction 

A molecular spectrum arising from electronic transitions 
consists of many lines that correspond to different initial 
and final vibrational and rotational states. Lines are also 
broadened and shifted by various mechanisms making 
details of the spectrum difficult to predict. However, for 
some purposes, the interest is primarily in the total radia- 
tion or in a rough idea of local features. In such cases, 
calculations can be considerably simplified. We will derive 
several approximate formulas connected with different 
models of the spectrum which are useful in problems of 
this kind. 

We begin by defining notation and reviewing some basic 
relationships. The energy levels of a molecule will be 
denoted as shown in Fig. 1. A given electronic state is 
composed of many vibrational states (with energy E v above 
the lowest vibrational level) and rotational states (with 
energy Ej above the lowest rotational level in a given 
vibrational state). As shown in Fig. 1, the lowest level in 
the electron state has energy E e above the ground state of 
the molecule. 


If the molecules are in thermal equilibrium and obey 
Boltzmann statistics, the fraction of all molecules in a par- 
ticular ev] level is 


/rtT i 1\ T (Ee + Ev + Ej)- 1 

o^-< 2 ,+ 1 > e n — st — J 

N Q 


where g e is the electron degeneracy given by g e = 2S + 1 
for 2 states and 2 (2S + 1) for all other states. Here, S is 
the spin quantum number of the electron state. The 
denominator 


+ (2 ) 


is the total partition function. 
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We now write the rotational energies as 



The fraction of molecules in a given electron state, 
regardless of the vibrational or rotational level, is 



2 N ev j 
__ v, j 

N 


Ej = J(J + l)hcB e (5) 

where B e is a constant depending only on the electron 
state, and not on the vibrational state. This approximation 
is accurate enough for our purposes (Ref. 1). The rota- 
tional state with the highest population corresponds to 

d /AF . T J (J + 1) hcB e ~] n 

rfJ< 2; + i )exl, L If J = 0 


or 



Supposing the gas to be hot enough that J max > > 1, the 
rotational partition function 

Qj= y j (2/ + l)exp(-|j^ (6) 

J 

can be approximated by an integral 

/■- r J(J + l)hcB e l , w kT 

Q J =J 0 ( 2J + ex P [ kT \ dJ ~ hcB e 


Since we have assumed the rotational constant to be inde- 
pendent of the vibrational level, the total partition func- 
tion can be factored, and Eq. (4) can be written as 




N, 


N 


= F e 


exp (“ (2/ + !) ex P ( - 


Qv 


( 8 ) 


( 3 ) 


where the vibrational partition function is 


Equation (1) can then be rewritten as 


N ev j 

IT 


= Fe, 


exp ( _ §) (2/ + ^ exp ( _ §) 

2 ***{ - §) X) m + l)e,p ( ■ §) 

v j 

(4) 


g ’ = 2 eX P ("§) (9) 

V 

Suppose a molecule makes a transition from an upper 
state n to a lower state m. If A nm is the Einstein transition 
probability, hc\r nm the transition energy (with the wave 
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number ~ nw in cm 1 ), and the number density of molecules 
in the upper state N n) the total emitted intensity is 

I™ = N n hcv nm A nm (10) 

The transition probability is related to the emission oscilla- 
tor strength /JJJ by 


which shows the physical meaning of The approxi- 
mate separation of the oscillator strength into electronic 
and vibrational parts is a conventional procedure, but is 
not always very accurate. The oscillator strength will be 
separated in the following formulas, but the two parts can 
be easily recombined if desired. 

Using Eqs. (12), (8), and (14), the following is obtained: 


A 


nm 


8tt 2 

me 


fnm 
/ em 


( 11 ) 


fe'v'J'->e”v"allJ" 

N 


— 4rr’2hc 2 irr 0 r 3 j, _ 


f elec 
em 


qw " 


thus, defining r„ = e 2 /rnc~ = classical electron radius. 


l nm = N n • 4w • 2 he 2 7rr„T 3 nm f™ (12) 



(2/' + l)exp(-§) 


(16) 


The intensity per steradian is obtained by dividing by 4 tt. 
If all quantities are expressed in cgs units, the right-hand 
side must be multiplied by 10~ 7 to convert to the common 
units of W/cm 3 . Numerically, the total radiated intensity 
of the m transition per cubic centimeter of gas can 
then be written 

I nm = L32 x 10 - 23 v 3 nm N n f™ W/cm 3 (13) 


According to the selection rule A J = 0, ±1, three (or some- 
times two) lines are involved. Therefore, v], must be an 
average for these lines arising from a given /' level. They 
are actually so close together that for our purposes no 
important ambiguity is raised. 

It is sometimes convenient to rewrite Eq. (16) in terms 
of properties of the lower electron state. From Eqs. (3), 
(6), and (9) we have 


It should be emphasized that in these formulas, and in 
those that follow, a negligible self-absorption is assumed 
so that only optically thin layers of radiating gas are 
considered. 


Kt = 


geQ'vQjex p 


(- 5 ) 


The oscillator strength can be written approximately 
(Ref. 1) as the product of an oscillator strength f e e l ^ c for the 
electron transition as a whole, and a Franck-Condon factor 
q V ’ V - between given vibrational levels. Therefore, 

f n eZ=ft l : C qv'v" (14) 


g'e'Q'v'Q'Jex? 


F" 

*e,r 


(-§) 


(17) 


Also, emission and absorption oscillator strengths are 
related by 


is related to the intensity of transitions from a completely 
specified upper state (e',t/,J') to lower states having 
definite electron and vibrational levels (e" and t/'), but all 
allowed rotational levels J". Franck-Condon factors obey 
the sum rule 


21 qvv" 1 


g * i felec — p ff felec 
« / em &e / abs 


(18) 


so one obtains (Fig. 2) 

felecpr — exp I” _ E e ~ E e 1 c 

Tern r e,T Q' V 'Q'/ I JtT J V e > T 

Note also from Fig. 2 that 


(19) 


thus, if one sums over all final vibrational levels, 

V f nm = felec V n — felec 
^ / em / em ^ Hv'v" — / em 

v " V" 


(15) 


hcv r = E' e + E' v + E'j - (£" + E" + E'/) 
or 

K - E7 + E' v + E'j = hcVj. . + E' v ' + E7 (20) 
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Therefore, Eq. (16) can be rewritten as 

Je'v'J' -»■ e"v" all J” 

N = 


( hc~j - -*j-\ 
kT — ) 

esp (~w) (21 ' + 1) (^w) 


X q v x y 


Q" 


(21) 


original total intensity in any frequency interval contain- 
ing many lines. For given t/ and u", we define v v 0 v " by 
hcv%' v " = E' + E'- E" — E". Since rotational states have 
energy J (J -f 1) hcB e above the ] = 0 state, it follows that 
the wave number difference between the /' J" line and 
is 

ir - vf *" = y (/' + 1) B; - J" (/" 4- 1) B" (22) 

Since J can change by at most one unit in the transition, 
and since we have previously assumed that large /-values 
predominate, we approximate /' ^ J" ss J, so that 

+ i)<b;-b") (23) 

Smoothing out the rotational structure means that we 
treat J as a continuous variable. Differentiation of Eq. (23) 
gives dv — (2/ + 1) (B' — B") dj. We then spread the J th 
line over the interval dj — 1 (i.e., just smoothing it out to 
meet the neighboring line J + 1), so that the range of wave 
numbers covered by the line is 

Av = (2/ + 1) | Be — Be' | (24) 

using absolute value signs to make Air positive, by con- 
vention. 


Now, let faArbe the intensity in the wave number inter- 
val Air. From Eqs. (23) and (24), we have 


(2 /+ l)exp^--^r) = 

(2/ + 1) exp f ^ J - 


\b: 


At- r . 

7 - g" J ex P[~ 




B' 


kT (B' e — B") 


(25) 


so that Eq. (21) becomes 


II. The Smoothed Rotational Model 

The first model we will discuss is the smoothed rota- 
tional model, in which individual lines are smeared out 
to obtain a continuous spectrum. That is, the set of lines 
belonging to a given vibrational transition will be spread 
out over a continuum of frequencies, while retaining the 


hrAv = 


, OL 2 3 / hcv\„ m? hcB " 

4,r ■ 2 he Ttr 0 v exp ^ jNF e T g,, fcr | ^ _ B » | 


EV + 


X exp I 


he B 

b: - b: 




kT 


Av 


(26) 
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At a given wave number, several vibrational transitions 
may contribute, so one must sum over v' and v"\ 


I&ai A? = 


4it • 2hc 2 irr< t v e exp ( - Jj) NF "r^7r 


Z 


X > < 7 „ exp 

1 

V'V " 


hcB” , 

E V+ K _ BV (*-<’ v ") 

kT 


BV | 

A'v 

(27) 


Summing over those values of v' and v" for which 
'C'" — V (if B' e > B") 

K' r — V (if B' e < BY) 

It is also common to treat the intensity as a function of 
wavelength. Using I\&\ — Iv A~,v= 1/A, and A7 = A A/A 2 , 
it follows that 


7f 0(a! AA = 


1 / hr \ f elee 


licit" 


Z 


X > exp 

f 

v'v" 


b;'| 

_hcK (1 1 v 

- B7\ A XX*-) 


kT 


summed over v ' and o" for which 
Ap*^A (if 

< v "^a (if b;<b") 

These formulas for the smoothed rotational model have 
been derived before in Refs. (2-4). There are differences 
in notation between the various results, and some actual 
differences that are not important in the temperature and 
wavelength range assumed here (v/T > 2). 

We have programmed Eq. (28), and part of the com- 
puted spectrum of the CN violet transition (B 2 X-» X 2 l) 
is reproduced (see Fig. 4). A numerical integration over 


the spectrum can then be performed. The total radiated 
intensity per molecule (at 7000°K) for this transition is 
shown in Table 1. 


III. The First Vibrational Model 

The smoothed rotational model just presented is highly 
accurate for calculating total radiation or the averaged 
local details of the spectrum. However, considerable labor 
is required to apply this method to particular cases— 
perhaps more than is justified by knowledge of the physical 
parameters or the desired accuracy of the result. Fortu- 
nately, it is possible to derive simpler formulas that cor- 
respond to less detailed spectral models. 

To obtain the total radiated intensity per molecule in 
going from electron state e ' to electron state e " (i.e., includ- 
ing all bands), one can sum Eq. (16) over t/, v", and /". 
The result is 

""--“-^EE* -K-5) 


N 


X 


2> + D«p (-#>. 


(29) 


Expressing in terms of densities in the lower electron state 
(Eq. (19) along with the definition E' — E" = 0 from 

Fig. 1) can be used to obtain an alternate formula 


JAA 

— / l ji “ ^Trr Iab * e > T r 

N 4tt 2 hcnr 0 

(28) 

X p 


X 


(-§) 

(2 /' + 1) exp ( - -j^r) 


V' V" 


(30) 


As a simplification, we now neglect all rotational struc- 
ture, and assume, that to a good approximation for given 
vibrational states, all light is emitted at a single average 
frequency. That is, let 

(2/' + 1) exp(- = 

J ' 

(2/' + l)exp (- 

(31) 
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where tW' is an average wave number between vibra- 
tional states t/ and v”. This procedure is equivalent to 
using the simplified spectral model shown in Fig. 3b, as 


compared with the more detailed smoothed rotational 
model illustrated in Fig. 3a (also see Fig. 4). Then, using 
Eq. (29), the intensity per molecule for the complete elec- 
tron transition becomes 


v'-v" 



(a) SMOOTHED ROTATIONAL MODEL 



telec J7' 

4n • 2ihc 2 irr o 7fm C,T 
Vv 



q V v exp 




(32) 


or alternatively, from Eq. (30) 


Je' -+ c" 


4 tt • 2hc 2 Trr f lec F ” t ®' J exp ( - ^2!L\ 

4tt Zhc Trr 0 j abs exp ^ j 

X qvv exp ^ ~ 3 V , V „ (33) 


0-0 



(b) FIRST VIBRATIONAL MODEL 


It is comparatively simple to use these formulas to calcu- 
late the total intensity; since it is not necessary to consider 
all /-values, as in Eqs. (29) and (30); nor is it necessary to 
perform a numerical or graphical integration, as in the 
smoothed rotational model. The result of Eq. (33) for the 
CN violet transition is given in Table 1. 

IV. The Second Vibrational Model 

As a further simplification, suppose that to sufficient 
accuracy we can group the lines from a given upper vibra- 
tional state v ' around some average wave number v v >. 
That is. 


0-0 

0-1 1-0 



(c) SECOND VIBRATIONAL MODEL 

Fig. 3. Schematic spectra for three molecular 
band models 


yy qvvK’v = </»■»■ - v%. (34) 

V " V" 

using the sum rule 

2 qv’v" i 

V" 

for the Franck-Condon factors. From Eqs. (32) and (34) 
one can then write 


Te'^e" p' 

-^-=4 T -2^W^-g7- 


X 


^ Vexp (- A ) (35) 
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or alternatively, using Eq. (33) 


VI. Discussion 


Je' -»e" 

N 


4TT'2hc 2 irr 0 fll‘ c F‘ 


abg L e,T 


& 

Q'/Q'S 


Xe * p (-Tsr)X/' exp (“'i r ) 

(36) 

These formulas are clearly less accurate than those of the 
previous section, since it is more difficult to choose an 
average wave number However, they provide a con- 
siderable simplification for practical calculations, since the 
Franck-Condon factors are not needed. Often these are, 
at best, only poorly known. Therefore, it may be useless 
to attempt a more accurate calculation than suggested by 
Eqs. (35) and (36). This second vibrational model corre- 
sponds to the spectral model illustrated in Fig. 3c. 


V. The Simplest Model 

As a final model, we neglect structure altogether, and 
assume that all light from the electron transition (includ- 
ing all bands) is emitted at a single average frequency. 
That is, let 


X K exp (~ §) = X exp ( - w) =’'*■■« 


(37) 


where ~ av is the average wave number. The intensity for- 
mulas of Eqs. (35) and (36) then become 





or alternatively 


(38) 


Je'^e" O' O' 

= 4tt • 2hc 2 7rrJ e lec F" -± J ^ V PY n 

N abs " e, T q’jq:: exp 


/ hcv 00 \ 

\ ~w~ p™ 


(39) 


This monochromatic approximation is essentially that 
introduced in Ref. 5. It is the simplest, but least accurate 
of the models, because it is difficult to choose the proper 
value of 7 ai? . If a spectrum is not too diffuse, the value can 
be guessed, or simply chosen to be v av = ~ 0 o as defined in 
Fig. 1. This value was used for the CN violet transition, 
and the result is given in Table 1. 


An advantage of the simpler models is that they can 
be quite readily calculated by hand. Besides, the poorly 
known oscillator strengths, the Franck-Condon factors, or 
the rotational constants may make more elaborate calcula- 
tions pointless. However, one would like to know how to 
estimate the accuracy of the various models. Therefore, 
they were all programmed to provide direct comparisons 
in specific cases. In particular, each method has been used 
to compute the total intensity per molecule of CN violet 
radiation, corresponding to the B 2 2 X 2 2 transition. The 

results are shown in Table 1. 


Table 1. Total band CN violet radiation at T = 7000° K 


l/N X 10“ W/molecule 

Simple 

model 

Second 

vibrational 

model 

First 

vibrational 

model 

Smoothed 

rotational 

model 

Exact 

model 

1.87 

2.18 

2.16 

2.19 

2.20 


The oscillator strength used was f e a l e g c ~ 0.027, as re- 
ported in Ref. 6. The Franck-Condon factors were taken 
from Ref. 7, and the other molecular parameters are from 
Ref. 1. In the simplest model as derived in Section V, the 
v av value was chosen to be ~ 0() , as defined in Fig. 1, rather 
than attempt to estimate an average wave number. Simi- 
larly, for the second vibrational model, v v r was chosen to 
be the wave number of the v —v ff transitions, and for the 
first vibrational model, ~ V ' V " was chosen to bet^' v '' as used 
in Eq. (22). 


The smoothed rotational model provides an actual 
detailed spectral distribution. The computer found the 
intensity as a function of wavelength, a sample of which 
is shown in Fig. 4. The spectrum was then integrated to 
obtain the total intensity. Finally, the total radiation was 
calculated by directly evaluating the sums in Eq. (30). This 
latter result is quoted in the last column of Table 1. It will 
be observed from the table that there is good agreement 
between the models, even for the simplest model, indicat- 
ing that it would have been adequate for many applica- 
tions to use one of the simpler models to find the total 
intensity. 


Figure 5 shows the spectral distribution obtained in the 
first vibrational model, and compares it with an exact 
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I/N , W/molecule-cm 



WAVELENGTH, A 

Fig. 4. CN violet bands calculated with the smoothed 
rotational model (individual v'-v" bands 
are identified) 


calculation using Eq. (21). The results of the latter com- 
putation were smeared out to give a continuous spectrum. 
As shown in the figure, the first vibrational model pro- 
vides only a rough idea of local features, but the model 
gives a good result for the total integrated intensity, as 
Table 1 indicates. 



WAVELENGTH, A 

Fig. 5. Comparison of two molecular band models 
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